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ABSTRACT 

The unipolar induction DC circuit model, originally developed by Goldreich & Lynden-Bell for the 
Jupiter-Io system, has been applied to different types of binary systems in recent years. We show 
that there exists an upper limit to the magnetic interaction torque and energy dissipation rate in 
such model. This arises because when the resistance of the circuit is too small, the large current flow 
severely twists the magnetic flux tube connecting the two binary components, leading to breakdown of 
the circuit. Applying this limit, we find that in coalescing neutron star binaries, magnetic interactions 
produce negligible correction to the phase evolution of the gravitational waveform, even for magnetar- 
like field strengths. However, energy dissipation in the binary magnetosphere may still give rise 
to electromagnetic radiation prior to the final merger. For ultra-compact white dwarf binaries, we 
find that DC circuit does not provide adequate energy dissipation to explain the observed X-ray 
luminosities of several sources. For exoplanetary systems containing close-in Jupiters or super-Earths, 
magnetic torque and dissipation are negligible, except possibly during the early T Tauri phase, when 
the stellar magnetic field is stronger than 10 3 G. 

Subject headings: binaries: close - stars: magnetic fields - stars: neutron - white dwarfs - gravitational 
waves - planetary systems 



1. INTRODUCTION 

In a seminal paper, Goldreich & Lynden-Bell (1969) 
developed a DC circuit model for the magnetic interac- 
tion between Jupiter and its satellite Io (see also Pid- 
dington & Drake 1968). In this model, the orbital mo- 
tion of Io in the rotating magnetosphere of Jupiter gen- 
erates electromotive force (EMF), driving a DC current 
between Io and Jupiter, with the closed magnetic field 
lines serving as wires. In essence, the Jupiter-Io system 
operates as a unipolar inductor. This model helped to 
explain Jupiter's decametric radio emissions which are 
correlated with the motion of Io and the hot spot in 
Jupiter's atmosphere that is linked magnetically to Io 
(e.g., Clarke et al. 1996). 

In recent years, the DC circuit model of Goldreich & 
Lynden-Bell has been adapted and applied to other types 
of astrophysical binary systems, including (i) Planets 
around magnetic white dwarfs (Li et al. 1998); (ii) Ultra- 
compact white dwarf binaries with periods ranging from 
a few minutes to a hour (Wu et al. 2002; Dall'Osso et 
al. 2006,2007); (hi) Coalescing neutron star - black Hole 
binaries and neutron star - neutron star binaries prior to 
the final merger (McWilliams & Levin 2011; Piro 2012; 
see also Vietri 1996, Hansen & Lyutikov 2001 and Lyu- 
tikov 2011 for more general discussion of unipolar induc- 
tion in such binaries); (iv) Exoplanetary systems con- 
taining short-period (~ days) planets orbiting magnetic 
MS stars (Laine & Lin 2012); 

In the DC circuit model, the magnetic interaction 
torque and the related energy dissipation rate are in- 
versely proportional to the total resistance T^tot of the 
circuit, including contributions from the two binary com- 
ponents and the magnetosphere. Thus it seems that ex- 
treme power may be produced for highly conducting bi- 
nary systems. However, we show in this paper that when 



T^tot is smaller than a critical value (~ 480v re i/c Ohms, 
where v re \ is the orbital velocity as seen in the rotating 
frame of the magnetic star), the large current flowing in 
the circuit will break the magnetic connection between 
the binary components. Thus there exists an upper limit 
to the magnetic torque and the energy dissipation rate. 
In Sects. 2-4, we examine the general aspect of the DC 
circuit model, the upper limit and its connection with 
"open circuit" effect ("Alfven radiation"). We then dis- 
cuss applications of the model to various astrophysical 
binaries in Sects. 5-8 and conclude in Sect. 9. 

2. DC CIRCUIT MODEL OF MAGNETIC INTERACTION 

Consider a binary system consisting of a magnetic 
star (the "primary", with mass M*, radius R+, spin ft s , 
and magnetic dipole moment /i) and a non-magnetic 
companion (mass M c , radius R c ). The orbital separa- 
tion is a, and the orbital angular frequency is ft. The 
magnetic field strength at the surface of the primary is 
B* = fJi/R+. The whole binary system is embedded in 
a tenuous plasma (magnetosphere). For simplicity, we 
assume fi, Cl s and /j, are all aligned. 

The motion of the non-magnetic companion relative to 
the magnetic field of the primary produces an EMF £ ~ 
2R C \E\, where E = v re i x B/c, with v re i = (ft — Q s )a(p 
and B = (—ja/a 3 )z. This gives 
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where Aft = ft — ft s . The EMF drives a current along 
the magnetic field lines in the magnetosphere, connecting 
the primary and the companion through two flux tubes 
(Fig. 1). The current in the circuit is 



1 = 



(2) 
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Fig. 1. — DC circuit model of magnetic interactions in binary- 
systems a la Goldreich Sz Lynden-Bell (1969). 

where the total resistance of the circuit is 

T^tot = + He + 27^ mag , (3) 

with 1Z C1 TZmzg the resistances of the magnetic star, 
the companion and the magnetosphere, respectively. The 
energy dissipation rate of the system is then 



diss 



21 2 R t , 
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tot 
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where the factor of 2 accounts for both the upper and 
lower sides of the circuit. 

Note that because of the finite resistances 72* and 1Z C1 
the magnetic flux tube in general rotates at a rate Qf 
between Q s and Q. The slippage of the companion rela- 
tive to the magnetic field gives £ c c± (2jiR c / co 2 ){Q j — Q j f)- 
The slippage of the primary's rotation relate to the flux 
tube generates ~ (2{iR c /ca 2 )(QF — Q s ) at the mag- 
netic foot on the surface of the primary. The total EMF 
(£* + £ c ) is the same as in Eq. ffl, independent of Op. 

The total magnetic force (in the azimuthal direction) 
on the companion is c± (2R c )(2TB z /c), with B z = 
— /i/a 3 . Thus the torque acting on the binary's orbital 
angular momentum is 



T = J rb — -aR c XB z 
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The torque on the primary is I*Cl s = — T (where I* is 
the moment of inertia). The orbital energy loss rate as- 
sociated with T is then 
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The total energy dissipation rate of the system is 
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in agreement with Eq. Q. 

3. MAXIMUM TORQUE AND DISSIPATION 

The equations in the previous section clearly show that 
the binary interaction torque and energy dissipation as- 
sociated with the DC circuit increase with decreasing to- 
tal resistance 7Z to t- Is there a problem for the DC model 
when T^tot is too small? The answer is yes. 



The current in the circuit produces toroidal magnetic 
field, which has the same magnitude but opposite di- 
rection above and below the equatorial plane. The 
toroidal field just above the companion star (in the up- 
per flux tube) is B<f,+ c± —(2it/c)J r , where the (height- 
integrated) surface current in the companion is JT^ 00 
— 4X/ (ttR c ). Thus the azimuthal twist of the flux tube is 
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where v re \ = aAQ = a(Q — O s ) is the orbital velocity in 
the corotating frame of the primary star. Clearly, when 
T^tot is less than 16v re i/c 2 , the flux tube will be highly 
twisted. 

GL already speculated in 1969 that the DC circuit 
would break down when the twist is too large. (For the 
Jupiter-Io system parameters adopted by GL, the twist 
IC0I <C 1.) Since then, numerous works have confirmed 
that this is indeed the case. Theoretical studies and nu- 
merical simulations, usually carried out in the contexts 
of solar flares and accretion disks, have shown that as a 
flux tube is twisted beyond <: 1, the magnetic pressure 
associated with B^ makes the flux tube expand outward 
and the magnetic fields open up, allowing the system to 
reach a lower energy state (e.g., Aly 1985; Aly & Kui- 
jpers 1990; van Ballegooijen 1994; Lynden-Bell & Boily 
1994; Lovelace et al. 1995; Uzdensky et al. 2002). Thus, 
a DC circuit with ^ 1 cannot be realized: The flux 
tube will break up, disconnecting the linkage between the 
two binary components. 

A binary system with 7£ tot <J Wv re \/c 2 cannot estab- 
lish a steady-state DC circuit. The electrodynamics is 
likely rather complex. At best we can expect a quasi- 
cyclic circuit, involving several steps: (a) The magnetic 
field from the primary penetrates part of the compan- 
ion, establishing magnetic linkage between the two stars; 

(b) The linked fields are twisted by differential rotation, 
generating toroidal field from the linked poloidal field; 

(c) As the toroidal magnetic field becomes comparable 
to the poloidal field, the fields inflate and the flux tube 
breaks, disrupting the magnetic linkage; (d) Reconnec- 
tion between the inflated field lines relaxes the shear and 
restore the linkage. The whole cycle repeats. 

In any case, we can use the dimensionless azimuthal 
twist to parameterize the magnetic torque and energy 
dissipation rate: 
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(10) 



The maximum torque and dissipation are obtained by 
setting ^ 1. Note that in the above, T is negative 
since we are assuming Q > Q s . A reasonable extension 
would let ((f) = C(Aft)/ft, with ( > 0. 

4. ALFVEN DRAG 

As discussed in GL, the validity of the DC circuit 
model requires that the slippage of the flux tube rela- 
tive to the companion during the round-trip Alfven travel 
time (tjC) along the flux tube be much less than R C1 i.e., 
(fi — VLp)atA <^ R c - When this condition is not satis- 
fied or when the poloidal field opens up, the disturbance 
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generated by the companion's orbital motion v re i will 
propagate along the field line as Alfven waves and radi- 
ate away (Drell, Foley & Ruderman 1965). The Alfven 
radiation power associated with the "open circuit" is 



1 



1 



. (B z v rel /v A yv A (27rR 2 c ) = -B 2 z R 2 c ^, (11) 
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where v A is the Alfven speed in the magnetosphere (as- 
suming v re \ ^ va)- The associated drag force is simply 
— —EA\i/vveh and the torque on the orbit is 
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~B 2 z R 2 c a V -^. 
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Comparing Ea\{ with Eq. (10), we find 
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Thus, ^Aif is always smaller than the maximum £?diss of 
a DC circuit. Equation ([To]) (with ~ 1) represents the 
maximum magnetic dissipation rate of the binary system, 
regardless of the details of the electrodynamics. 

5. NEUTRON STAR - NEUTRON STAR BINARIES 

Gravitational wave (GW) emission drives the orbital 
decay of the neutron star (NS) binary, with timescale 
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where in last equality we have adopted = 1.4M 
and mass ratio q = M c /M* = 1. The magnetic torque 
tends to spin up the primary when Q > Q s . Spin-orbit 
synchronization is possible only if the synchronization 
time t syn = 7*0/ |T| is less than tow at some orbital 
radii. With I* = kM*R 2 and the torque ([9|, we find 
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where in the second line we have adopted k = 0.4 and 
R* = R c = 10 km. Clearly, even with magnetar-like 
field strength {B* ~ 10 15 G) and maximum efficiency 
(((f) ~ 1), spin-orbit synchronization cannot be achieved 
by magnetic torque. [It was already known that tidal 
torque, due to both equilibrium tide (Bildsten & Cutler 
1992) and resonant tides (Lai 1994), cannot synchronize 
the NS spin during binary inspiral.] 

For the same reason, the effect of magnetic torque on 
the number of GW cycles during binary inspiral, iV, is 
small. We find 
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where / = Q/tt is the GW frequency, and 
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is the usual leading-order point-mass GW cycles. The 
correction factor due to the magnetic torque is 
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where /i m = M+q/(l + q) is the reduced mass of the 
binary. With Eqs. ([□])- (15]), we see that GW phase 
error a(dN/d\n /)q is much less than unity even for 
B* ~ 10 15 G and maximum ~ 1. 
The energy dissipation rate is 
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where in the second line we have used v re \ ~ rQ (for 
Q s <C fl) and adopted canonical parameters (M* = M c = 
1.4M , R+ = R c = 10 km). The total energy dissipation 
per lna is 
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Some fraction of this dissipation will emerge as elec- 
tromagnetic radiation counterpart of binary inspiral. 
Whether it is detectable at extragalactic distance de- 
pends on the microphysics in the magnetosphere, includ- 
ing particle acceleration and radiation mechanism (e.g., 
Vietri 1996; Hansen & Lyutikov 2001). 

Piro (2012) recently applied the DC circuit model to 
NS binaries. In the case when the magnetic NS primary 
dominates the total resistance, he found that the mag- 
netic torque synchronizes the NS spin prior to merger and 
significantly affects the gravitational wave cycles, even 
for modest (~ 10 12 G) NS magnetic fields. Using the 
resistance computed by Piro, we find from Eq. ([sj) that 
the corresponding azimuthal twist is much larger (by 
a factor ~ 10 8 at Q — 100 s _1 ) than unity, violating the 
upper limit discussed in Sect. 3. 

If one assumes that the magnetosphere resistance is 
given by the impedance of free space, 7£ mag = 47r/c, then 
the corresponding twist is = 2v re i/(7rc), which satisfies 
our upper limit. The energy dissipation rate is then 
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Not surprisingly, this is approximately the same as the 
Alfven power Eau [ see Eq. ([13|] with va = c and is in 
agreement with the estimate 01 Lyutikov (2011). 

6. NEUTRON STAR - BLACK HOLE BINARIES 

The situation is similar to the case of NS-NS bina- 
ries. In the membrane paradigm (Thorne et al. 1986), a 
black hole (BH) of mass Mh resembles a sphere of ra- 
dius R c = Rh = 2GMh/c 2 (neglecting BH spin) and 
impedance 1Zh = 47r/c. Neglecting the resistances of the 
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magnetosphere and the NS, the azimuthal twist of the 
flux tube in the DC circuit is 



(22) 



which satisfies our upper limit (Sect. 3). The energy 
dissipation rate is (cf. Lyutikov 2011; McWilliams & 
Levin 2011) 
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10 13 Gy \10mJ \3R h 
where we have assumed Mbh/M* ^> 1. 
7. ULTRA-COMPACT DOUBLE WHITE DWARF BINARIES 

Wu et al. (2002) and DalPOsso et al. (2006,2007) (see 
also Wu 2009) developed the DC circuit model for ultra- 
compact double white dwarf (WD) binaries, particularly 
for the systems RX J1914+24 (period P = 9.5 mins) 
and RX J0806+15 (P = 5.4 mins). Usual mass trans- 
fer models appear to have difficulties explaining some of 
the properties of these systems (e.g., the observed or- 
bital decay). The DC circuit model seeks to account for 
the observed X-ray luminosity (10 35 — 10 36 erg s _1 for 
RX J1914+24 assuming a distance of 100 pc) without 
mass accretion, while allowing for orbital decay driven 
by gravitational radiation. 

Using Eq. (10) with parameters appropriate to com- 



pact WD binaries, we find 
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Note that \i = 10 32 Gem 3 corresponds to B* ~ 0.5 MG 
(for = 6000 km), approaching the field strengths of 
Intermediate Polars. Obviously, even with the maximum 
asynchronization (AQ/Q = 1) and maximum efficiency 
((cf) ~ 1), i^diss falls far short of the observed X-ray lu- 
minosities. Wu et al. (2002) calculated the resistance of 
the WD and used Eq. Q to obtain a much higher energy 
dissipation power (see their Fig. 3) - evidently, their re- 
sult (which was also adopted by DalPOsso et al. 2006) 
corresponds to ^> 1, violating our upper limit. 

8. CLOSE-IN EXOPLANETARY SYSTEMS 

Laine & Lin (2012) recently applied the DC circuit 
model to study the interaction of close-in planets with 
the magnetosphere of host stars. From their calculation 
of the resistances of the planet, host star and magne- 
tosphere, they suggested that magnetic interaction may 
affect the orbital evolution of close-in super-Earths on a 
few Myr timescale, and produce hot spots on the surface 
of the host stars. 

Applying Eq. pi to a planetary system (M c = M p: 
R c = R p ), we fina that the magnetic torque induces or- 
bital decay (assuming Q > as is the case for close-in 
planets in a few day orbit) at the rate 
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The timescale is 
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where p p is the mean density of the planet. Thus it is 
clear that even at maximum efficiency ((^> ~ 1), a hot 
Jupiter (R p ~ Rj) or super-Earth (R p ~ 0.1i?j) in a 
P ~ 3 d orbit (a ~ 0.04 AU) around a solar- type star 
(with a typical dipole field of a few Gauss) will suf- 
fer negligible orbital decay due to the magnetic torque. 
If super-Earths migrate to their close-in locations in the 
early, T-Tauri phase of the star, when the stellar mag- 
netic field is stronger (a few kG), they may experience 
appreciable orbital evolution due to magnetic interac- 
tions. 



The energy dissipation rate is [Eq. (10)] 
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If this energy is accumulated (e.g., building up twist from 
C</> ~ to (cf, ^ l) over time At ~ 2tt / Aft ~ 3 d and 
then released suddenly (~ hours), the energy release is ~ 
E di SS At/ 2 ~ 10 26 erg (for the canonical parameter values 
adopted in the above equation). This is much smaller 
than the energy release of solar flares (10 29 — 10 32 ergs) 
and of superflares from solar- type stars (^ 10 33 ergs; see 
Maehara et al. 2012). 

Finally, it is instructive to compare the magnetic 
torque T mag = T [Eq. with the tidal torque (due 
to tide raised on the star oy the planet). Parameterizing 
tidal dissipation by the quality factor Q*, we have 
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where = 3Q*/(2&2) and ki is the Love number (Gol- 
dreich & Soter 1966). Thus, 
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With ~ 10 6 — 10 8 , the magnetic torque generally 
cannot compete with the tidal torque. 

9. CONCLUSION 

Closed DC circuit connecting a magnetic primary and 
a nonmagnetic companion (Sect. 2) can be an efficient 
unipolar engine in a binary system, potentially more 
efficient than an "open" circuit engine (see Sect. 4). 
The power of this DC engine is inversely proportional 
to the total resistance 7£ t ot of the circuit. However, 
we have shown that when T^tot is less than a critical 
value, ~ 16v re \/c 2 (in the cgs units) [see Eq. ^] or 
480(v re i/c) Ohms, the magnetic flux tube connecting the 
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two binary components will be highly distorted and the 
circuit will break. In this case, a quasi-cyclic unipolar 
engine may operate in the system (Sect. 3). Thus, there 
exists an upper limit to the magnetic interaction torque 
and the associated energy dissipation rate [Eqs. ([9])- (p!Q|)] 
of any unipolar engine in magnetic binaries. Several pre- 
vious applications of the DC circuit model to different 
types of binary systems apparently violated this upper 
limit. We have shown that: 

(i) In coalescing double neutron star (NS) or NS - 
black hole binaries: Magnetic interactions cannot syn- 
chronize the NS spin with the orbital motion and have 
a negligible effect on phase evolution of the gravitational 
waveform, even for magnetar field strengths. Neverthe- 
less, the energy dissipation associated with the unipolar 
engine can produce electromagnetic radiation prior to 
binary merger, although the detectability of this radia- 
tion depends on the microphysics of the binary magne- 
tosphere. 



(ii) In Ultra-compact white dwarf binaries: Magnetic 
energy dissipation is too small to account for the ob- 
served X-ray luminosities. Thus, the puzzling behaviors 
of several sources (RX J1914+24 and RX J0806+15) can- 
not be explained by the unipolar inductor circuit model. 

(hi) In close-in exoplanetary systems: Interaction be- 
tween hot Jupiters or Super-Earths and the magneto- 
sphere of their host stars does not lead to appreciable 
orbital evolution, with the possible exception of the early 
T Tauri phase, when the stellar dipole magnetic field is 
higher than 10 3 Gauss. Magnetic energy dissipation in- 
duced by the orbital motion of planets is generally negli- 
gible compared to the observed energy releases in stellar 
flares or superflares. 
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